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I. INTRODUCTION 

The interactions between different quasiparticles in 
the solid often lead to interesting effects. For exam- 
ple, exciton-polaritons or cavity polaritons appear in 
the microcavitiesi as the elementary excitations formed 
by the coupling of excitons and photons. Having ex- 
tremely small effective mass (10~ 4 — 10~ 5 of the effec- 
tive mass of the free electrons), obeying bosonic statis- 
tics and efficiently interacting with each other, polaritons 
can undergo a transition into a collective phase coherent 
state usually referred as polariton BEC2 accompanied by 
the onset of superfluidity 3 -. Some indications of the oc- 
curence of this effect in real systems have been recently 
reported^—. Another candidate for the achievement of 
the BEC in the domain of the condensed matter is a sys- 
tem of cold indirect excitons^. 

Another type of elementary excitations in the solid are 
phonons - the collective oscillations of the crystal lat- 
tice that can also interact with light. At relatively small 
light intensity emmition or absorption of phonons results 
in Mandclstam-Brillouin scattering with spectral energy 
shift on the phonon energy. However with the increase of 
the intensity of the light the probability of the scattering 
process grows and it can become comparable with the 
lifetime of the scattered states. As it was shown in^r— 
the strong reconstruction of the bulk phonon-polariton 
spectra occurs and the dispersions of both polaritons and 
phonons appear to depend strongly on the light intensity. 

The interactions between the exciton-polaritons and 
strong acoustic waves has been a subject of recent 
studies^ showing several interesting effects and open- 
ing many promising directions in the future, such as the 



possibility to observe polariton Bloch oscillations 13 . 

In the present paper we consider the interaction be- 
tween a condensate of quantum well (QW) excitons or 
cavity exciton-polaritons and a coherent phonon field, 
possessing the same dimensionality (2D or ID). The co- 
herent 2D acoustic wave can be obtained either by an 
external pumping (as ini£) or by quantizing the phonon 
modes within an acoustic waveguide^ around the QW 
containing excitons (for cavity polaritons, the waveg- 
uide should be embedded inside the microcavity). If 
the phonons are not confined to the same dimensionality, 
they play a role of an incoherent reservoir providing ex- 
citon or exciton-polariton relaxation^, which is not the 
case we would like to consider here in order to obtain 
strong coupling. 

In the phonon-photon or phonon-exciton coherent in- 
teraction models developed before^, the phonon field 
has usually been treated as an external classical field. 
This approach is valid under strong external pumping, 
whereas the model we develop can be used both for the 
external pumping case and for the case when phonons 
are created by the excitons themselves. 



II. FORMALISM 

We treat our system consisting of three interacting 
classical fields: complex fields ip± corresponding to right 
and left circular polarized excitons or exciton-polaritons 
and a real vector field of the lattice displacements u corre- 
sponding to phonons. The energy of polariton or exciton 
system is measured from the minima of the dispersion 
curve. The Lagrangian of the system reads: 
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The first line corresponds to free polaritons described by 
means of their macroscopic wavefunctions ib s and acous- 
tic phonons described interms of the lattice displacement 
field u = {u X i u y ). The second line to their mutual inter- 
actions: the term with a\ describes interactions of the 
polaritons of the same circular polarization, term with 
a-2- interactions of the polaritons of the opposite circu- 
lar polarization, term with g- polariton-coherent phonon 
interaction. For the moment, we neglect the wavevector 
dependence of this interaction and g should be considered 
as a constant determined by the deformation potential. 
m is an effective mass of cavity polaritons and p and Y 
are the density and the Young's modulus, determining 
the sound velocity c s = \/Y/p. 

The equations of motion for the considered fields are 
Lagrange equations (ib s and ip* should be considered as 
independent functions): 



where, in order to have the most general description, we 
have introduced the terms corresponding to the resonant 
pumping of the polariton and phonon modes (P^ ans 
Pcf, respectively) and the finite lifetime of the polaritons 
and phonons described by terms containing 7^ and 7^ 
respectively. 

III. ANALYTICAL SOLUTION 

In this section we calculate analytically the dispersion 
of excitations of coupled exciton or exciton-polariton and 
phonon modes with a k-independent interaction g. For 
simplicity we neglect here the spin of the polaritons, as 
our main goal here is to investigate the effects of spin- 
independent polariton- phonon coherent coupling. Then, 
one has the system of the equations for two coupled fields: 
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Taking the divergence of the second equation and intro- 
ducing a scalar phonon field 



4> = div(u) 



(6) 



one gets 



z dib s 
dt 



h 2 
2m 



dib h 2 
ih^- = -i^ - — VV + a\il>\ 2 i> + gW + P(r, t) (9) 
dt 2m 



o 2 



dt 2 



V 2 [Ycb + gW 



Pl<t> d( t> 
h dt 



+ P ph (r,f)(10) 



Now, consider a spatially homogeneous resonant pump 
of the polaritonic field under normal incidence, P = 
P e~ iuot ,P ph = 0. Then, introducing ib = ^ e ~ iuiat one 
gets 
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Looking for a spatially homogeneous stationary solution 
*o,</>o one gets 



- (lo + ili,) *o + a|* | 2 * + .g0o*o + P = (13) 

Y^+g\^ \ 2 =F (14) 

where F is hydrostatic pressure in the absence of the 
polaritons. In what follows we assume it to be zero. 



0o = -9|*o| 2 /^ (15) 
and for determination of the polariton field one has: 
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- (w + *7^) *o + U - |*o| 2 *o + P = (16) 

This is the same equation as for the case with 
polariton-polariton interactions only but with a renor- 
malized polariton-polariton interaction constant, 



a = a- — (17) 

One sees, that as correction to the interaction constant 
is negative, polariton-polariton interactions mediated by 
acoustic phonons are attractive. This is not surprising, 
as similar interactions for the electrons give attractive 
BCS potential. 

The equation 1161 is well studied in context of quantum 
microcavities and can lead e.g. to the S-shaped depen- 



dence of the intensity of the polariton field on pump in- 
tensity, which results in a bistability and hysteresi a 18 ' 19 . 

Let us now calculate the dispersion of the weak el- 
ementary excitations. Following a standard procedure, 
one represents the solutions in a following form 

* = * + Ae 1 ^-^ + s*e-^ kr -"**) (18) 
<t> = 0o + Ce 1 ^-^ + c* e -*(kr-^t) (19) 

which should be then put into Eqs. 1121121 Consider- 
ing deviations from equilibrium being small, one should 
then carry out the procedure of linearization, i.e. in the 
resulting equations keep only the terms linear in A,B,C. 

The system of equations for determination of the small 
amplitudes A,B,C reads: 
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Thus the dispersions of the elementary excitations can be determined from the following equation: 
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Here we used the relation [15] between the amplitudes 
of the excitonic and phonon fields. 

The consideration of the incoherent continous pump 
of the polariton system accompanied by the onset of the 
polariton BEC should be made in a different way. In 
this case, one can assume that in a stationary regime the 
finite lifetime of the polaritons is compensated by the 
pump, and calculations can be done in the same way as 
for particles with infinite lifetime with no external coher- 
ent pumping (i.e. one puts P — 0,7^ = 7^, = 00). Be- 
sides, the energy of the macroscopically occupied mode 
is not pinned by the frequency of the resonant pump 
wo but is given by the chemical potential of the system 

p = [a — Sy^j l^ol 2 - In this case, the dispersions of the 
elementary excitations read 



where 8 = an + . 

This gives a following equation for u(k), the disper- 
sion of the new quasiparticles in the strongly-coupled 
system of bogolons (BEC elementary excitations )2£ and 
phonons: 

- (h 2 p)uj 4 + [pE (k) (E {k) + 2cm) + Yh 2 k 2 ] uj 2 +(23) 
+k 2 E (k) [2g 2 n - Y (E (k) + 2cm)] = 

with E (k) = 4^-. This equation gives rise to two solu- 
tions, close to the original dispersions in the low-fc limit 
and showing the typical anticrossing at higher wavevec- 
tors. As expected, in the region of small k the dispersions 
are sound-like, 
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FIG. 1: (Color online) Dispersion of elementary excitations 
of an exciton condensate strongly coupled with the phonon 
field. Dashed lines show the original excitations, solid lines 
show the renormalized dispersions. Inset: same dispertions 
for the small k-vectors 
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Figure 1 shows the dispersions calculated with the use 
of this equation (solid lines) in comparison with the orig- 
inal dispersions at zero exciton-phonon coupling g = 0. 
The inset demonstrates the linear character of the dis- 
persions near zero. Here and below we consider only 
exciton condensates, although the theory developed ap- 
plies to both exciton and exciton-polaritons (the latter in 
the parabolic approximation). The material parameters 
used are those of GaAs, as i n 22 ' 23 . Since the dispersion 
of phonons is not steep at all compared to that of polari- 
tons, the anticrossing of the two dispersions takes place 
at relatively small wave vectors (10 3 m _1 ). For cxcitons 
the situation is much more favorable, because their en- 
ergy grows slower with k than that of polaritons, and 
thus the anticrossing takes place at much larger k (10 7 
m -1 ). In order to demonstrate the effects linked with the 
presence of the condensate we have to choose its density 
carefu lly, so t hat the speed of sound in the condensate 
c = sj a\njm lies below that of phonons, otherwise the 
branches will not anticross. 



IV. WAVEVECTOR DEPENDENCE 

In order to obtain an analytical solution for the dis- 
persions of excitations, we had to assume that the in- 
teractions between excitons and phonons were indepen- 
dent on the wavevector. However, in realistic systems 
this is not the case: g(k) first increases with wavevector 
and then drops down to zero exponentially due to the 
overlap integrals between the exciton wavefunction and 
the amplitude of the acoustic wave. We will calculate 
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FIG. 2: (Color online) Dispersion of the excitations in the 
exciton condensate coupled with waveguided phonon mode. 
The wavevector dependence of interactions is taken into ac- 
count. Solid lines correspond to the stable situation with an 
additional valley appearing in the dispersion, whereas dashed 
lines show the unstable case with imaginary dispersion. 



this dependence for the structure placed in an acoustic 
waveguide. 

The matrix element of exciton-acoustic phonon inter- 
action is^: 



M ac (q) = 



hq 



2pc s SL 



G(q\\,<lz) 



(26) 



where S, L - dimensions of the QW, q - acoustic wave 
vector. One can represent the function G{q\\,q z ) as: 



G(q lb q z ) = Dj\\{%)li{<lz) +D h lH(q ll )Ijt(q z ) (27) 

D e , Dh - are the deformation coefficients for electrons 
and holes and I^{q z ), ^(^l); I h(%) are the 

overlap integrals between the phonon and exciton mode. 

In order to calculate these integrals we have to obtain 
the distribution of the amplitude of the guided phonon 
wave inside the waveguide. Taking into account that QW 
thickness is much less than the thickness of waveguide, 
we can assume that the amplitude of the acoustic field is 
constant over the QW in the z direction. This value can 
be found from the waveguide equations 2 ^: 
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where d is a thickness of the waveguide, n\, n% are 
acoustic diffraction coefficients of waveguide layers, and k 
is a parameter which can be found by solving an equation: 



k tan nd 



(29) 
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Putting the value of 4>q into the overlap integrals, we can obtain the final value of g(q\\). 
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First, we have calculated the dispersion of the excitations 
numerically, in order to compare the result with the an- 
alytical solution found above. The results of these cal- 
culations are shown in figure 2 (solid lines). One can 
see that if the maximum of the /c-dcpcndent exciton- 
phonon interactions coincides with the crossing of the 
original branches, the interactions can lead to a much 
larger splitting and even to the appearance of a valley 
on the dispersion, similar to the "roton minimum" 25 ' 26 . 
Dashed lines correspond to the extreme case, when the 
dispersion of the excitations becomes imaginary and the 
system becomes unstable against any small perturbation 
in this wavevector region. 

We have next studied the behavior of the system 
characterized by such interactions between bogolons and 
phonons. An anticrossing between two quantum levels 
can give rise to so-called Rabi oscillations, which take 
place if only one of the original levels is excited initially. 
We demonstrate them in our system by solving numeri- 
cally the equations for the exciton condensate order pa- 
rameter and for the displacement field (|12j). assuming 
a homogeneous condensate and a propagating acousting 
wave in the x direction as the initial conditions. The 
system is homogeneous in the y direction and thus the 
problem is reduced to ID. 
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FIG. 3: (Color online) Rabi oscillations between the excita- 
tions of the exciton condensate (bogolons) and phonons: the 
amplitude of the propagating waves in each component oscil- 
lates in time, as expected. Panel (a) shows the condensate 
density as a function of coordinate and time, panel (b) shows 
the amplitude of the acoustic wave. 

Figure 3 obtained this way demonstrates the Rabi os- 
cillations between the two types of propagating waves: 
the bogolons (panel a) and the phonons (panel b). Ini- 
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FIG. 4: (Color online) Exponential growth of a weak excita- 
tion in the unstable case. Panel (a) shows the evolution of the 
condensate density over time, panel (b) shows the amplitude 
of the acoustic wave over time. 



tially, the exciton condensate is homogeneous (no bo- 
golons), whereas the phonon field contains a single 
monochromatic propagating wave. This situation corre- 
sponds to an excited atom in an empty cavity in original 
Rabi oscillations. After 150 ps, the amplitude of the os- 
cillations of the acoustic wave drops to zero, whereas the 
condensate density exhibits a running wave. This cor- 
responds to the atom being in the ground state and a 
photon in the cavity mode, if one continues the analogy 
with Rabi oscillations. After about 300 ps, the acoustic 
wave becomes strong again while the bogolons disappear 
completely, which corresponds to one period of Rabi os- 
cillations. 

Finally, we have simulated the unstable situation cor- 
responding to the dashed lines in figure 3. An excitation 
is created in the phonon field with its wavevector in the 
region with imaginary dispersion of excitations (flat real 
part). Figure 4 shows the results of our simulations: the 
small initial perturbation (a monochromatic wave) grows 
over time, but does not propagate in space, as expected 
from the flat real part and positive imaginary part of 
the dispersion (see fig. 3). The perturbations grow in 
both components (condensate and phonon field, panels 
(a) and (b)). This interesting effect cannot be described 
with the previous models (such as reffl6h . because they 
do not allow to take into account the effect produced by 
excitons or photons on the acoustic field itself. 
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V. CONCLUSIONS 

We have studied the interaction of QW excitons and 
cavity cxciton-polaritons (condensed or uncondensed) 
with a coherent phonon field (guided acoustic mode). We 
have found an analytical solution and studied the renor- 
malized dispersion of excitations in different situations, 
demonstrating the anticrossing of the two branches, an 
appearance of a valley at k ^ 0. If the lower dispersion 



branch touches 0, the system becomes unstable against 
any small perturbation in the corresponding wavevec- 
tor range. We have performed numerical simulations to 
demonstrate the Rabi oscillations between bogolons and 
phonons in the stable regime and the exponential growth 
of the perturbations in the unstable regime. 
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(237252) and IRSES "POLAPHEN" (246912) projects. 
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